Velocity plateaus and jumps in carbon 
nanotube sliding 



Xiao-Hua Zhang", Ugo Tartaglino^ *, 
Giuseppe E. Santoro^'^, and Erio Tosatti^'"^ 

" Surface Physics Laboratory and Department of Physics, Fudan University, 

Shanghai 200433, China 

International School for Advanced Studies (SISSA), and INFM-CNR Democritos 
National Simulation Center, Via Beirut 2-4, 1-34014 Trieste, Italy 

International Centre for Theoretical Physics (ICTP), P.O.Box 586, 1-34014 

Trieste, Italy 



Abstract 

The friction between concentric carbon nanotubes sliding one inside the other has 
been widely studied and simulated, but not so far using external force as the driving 
variable. Our molecular dynamics (MD) simulations show that as the pulling force 
grows, the sliding velocity increases by jumps and plateaus rather than continuously 
as expected. Dramatic friction peaks (similar to that recently noted by Tangney et 
al. in Phys. Rev. Lett. 97 (2006) 195901) which develop around some preferential 
sliding velocities, are at the origin of this phenomenon. The (stable) rising edge 
of the peak produces a velocity plateau; the (unstable) dropping edge produces a 
jump to the nearest stable branch. The outcome is reminiscent of conduction in 
ionized gases, the plateau corresponding to a current stabilization against voltage 
variations, the jump corresponding to a discharge or breakdown. 



1 Introduction 



Several years ago Zettl and others showed that nanotubes can be made to 
slide one inside another [Tp|3|4] . The sliding friction proved to be rather low, 
and the inner tube could even execute weakly damped oscillations out and in 
the outer one. 
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Assuming nanotube sliding friction to imply no wear, the remaining sources 
of friction are through i) dissipation of phonons from and through the outer 
tube to the substrate onto which it is deposited, and ii) electronic friction 
processes, corresponding to any extra internal resistivity arising inside each 
tube due to scattering caused by the presence of the other tube [5]. 

Little is known about nanotube electronic friction, and only indirect evidence 
of its existence has been identified so far Phononic nanotube friction, which 
on the other hand is the focus of our attention here, has been theoretically 
studied by various authors, particularly using classical molecular dynamics 
(MD) simulations [TIMTO] . 

The main source of tube-tube phononic friction was so far believed to be the 
inner nanotube terminations |10j, while the "bulk" friction, arising from mu- 
tual sliding of interior tube walls, is often considered to be irrelevant. Because 
of that, many simulation studies did not attempt to describe realistically the 
wall-wall interactions, sometimes assumed of a crude Lennard- Jones (LJ) type, 
with parameters which may underestimate the tube-tube corrugation and thus 
the bulk friction. 

Recently Guo [8] remedied that aspect and focused attention on bulk nanotube 
friction, adopting more accurate interactions, such as those parametrized on 
ab-initio energies by Kolmogorov and Crespi [TT] . 

In a recent exciting development, Tangney et al. [12j showed that tube-tube 
sliding friction can grow very large for some special sliding velocity, for example 
V ~ llOOm/s. This finding was rationalized by invoking a resonance condition 
between the "washboard frequency" uj^ = ^irv/a (where a is the inter-tube 
potential corrugation period) and particular phonon frequencies of the outer 
tube (radial breathing and longitudinal acoustic) which happen to possess a 
group velocity approximately equal to the sliding velocity v. 

Methodologically, it would moreover be more interesting to describe nanotube 
friction as a function of the applied sliding force F, rather than of some ex- 
ternally imposed velocity v. Force is in fact the external control parameter 
in friction, whereas sliding velocity is rather an outcome than an input. Thus 
for example, phenomena like stick slip occur in the real situation of an exter- 
nally applied force, but they would not for a hypothetical externally applied 
velocity. 

We describe here the simulated behavior of nanotube (phononic) sliding fric- 
tion, now as a function of applied force F, pulling the inner tube against the 
outer one. We use realistic tube-tube interactions in otherwise standard clas- 
sical MD simulations. Because we aim at describing purely bulk friction, we 
slide infinite tubes one inside the other, without any terminations whatsoever. 
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Preliminary results of our simulations show several novel features. The inner 
tube sliding velocity does not grow smoothly for increasing pulling force. There 
are instead large plateaus where the velocity v stays nearly constant in spite of 
a growing force F. Plateaus terminate with velocity jumps, typically to another 
velocity plateau. A second aspect concerns rotation. Leaving the inner tube 
free to rotate, there is no proper rotation in nonchiral tubes, where there may 
still be pseudorotations. On the other hand, there is generally a true rotation 
of the inner tube in the chiral case. Here we find that the inner tube angular 
velocity also grows by plateaus and jumps, which occur in coincidence with 
those of the linear velocity. 

The plateaus and jumps are rationalized as the natural outcome of frictional 
peaks for growing velocity, the stable (rising) side of the peak leading to a 
plateau, the dropping (unstable) side leading to a jump, as depicted schemat- 
ically in Fig. [TJ The presence, in a general tube force-velocity characteris- 
tics, of several peaks accounts for the multiplicity of plateaus and jumps. For 
each given pair combination of sliding tubes, their frictional peaks need to be 
mapped out. The peaks are ostensibly due to parametric excitation of vibra- 
tional modes in the double wall nanotube caused by the inner tube sliding, 
but the details of these parametric excitation need to be worked out in more 
detail. 



2 Infinite nanotube sliding: Molecular dynamics simulations meth- 
ods 

The simulations presented here are restricted to two concentric nanotubes, 
a strongly commensurate case — where the inner tube is a (5,5) armchair, 
and the outer one is (10,10) — , and a nearly incommensurate case, with a 
chiral (11,2) nanotube sliding inside a (12,12) one, as illustrated in Fig. [2j 
The (5,5)-(10,10) simulation cell comprises 160 (inner) + 320 (outer) carbon 
atoms. The (11,2)-(12,12) simulation cell comprises 196 (inner) + 336 (outer) 
carbon atoms. The two concentric nanotubes have a common 2;-axis, but have 
no terminations: they are connected across the boundary of the simulation cell 
through periodic boundary conditions (PBC) in the z-direction. This implies 
that we concentrate on bulk sliding, which is naturally the relevant source of 
the frictional peaks [12]. 

The inner tube is pulled along the axis with a force F acting on each atom, 
while the outer tube is subject to two constraints: one that fixes the position 
of its center of mass, and another that prevents the rotation around the z-axis. 
(These are global constraints: every single atom of the nanotubes is free to 
move.) A Berendsen thermostat [13j acts on the atoms of the outer tube alone, 
maintaining it at a temperature T, thus removing the frictional energy. This 
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is intended to mimic the behaviour of the real system, where the inner tube 
dissipates only through its interaction with the outer tube, in turn thermalized 
by direct contact with the substrate. 

The potential used in modeling intra-tube forces is a standard Tersoff-Brenner 
|14j . while inter-tube forces are described by the most recent Kolmogorov- 
Crespi potential [TT]. This choice was made following Ref. [8], who argued 
that a LJ interaction would underestimate the corrugation of the inter-tube 
potential, as is known to be the case in graphite [S]. 



3 Results 

Fig. [3] shows the time-dependence of the (5,5) inner-tube velocity when sliding 
inside a commensurate (10,10) outer tube under a fixed force F . The starting 
point was from rest. The outer tube is maintained at T = 300 K and impeded 
to rotate or translate globally. Two velocity plateaus are seen at f ?a 450 m/s 
and V ~ 780 m/s, separated by large velocity jumps. 

A similar occurrence of plateaus and jumps is observed for the incommensurate 
(11,2)-(12,12) case, shown in Fig. H] (here the outer tube is at T = 50 K), with 
the extra ingredient that the translational velocity plateaus and jumps are 
accompanied by a corresponding set of rotational plateaus and jumps in the 
global angular velocity of the inner (chiral) tube. 

We have just two examples here; but all indications suggest that these remark- 
able plateau-jump frictional phenomena are ubiquitous in nanotubes. In order 
to understand them we will concentrate here on the simpler commensurate 
(5,5)-(10,10) case, leaving the more complex chiral case for future studies. We 
carried out additional simulations, now at at fixed speed, to explore in more 
details the v-F characteristics. The simulations are performed by imposing a 
time-dependent force F{t), equal for all atoms, which enforces the constraint 
that the inner tube has to move with a given velocity f , following the general 
procedure of Ref. [T3] . 

After an initial transient, the force F{t) effectively balances (on average) the 
friction forces due to the outer tube, and provides, therefore, a direct mea- 
surement of -Ffriction, -^friction = F{t), the bar indicating a time-average. Fig. [5] 
shows the results obtained for Ffriction versus f , at temperature T=50 K. There 
is a clear peak of -Ffriction at f ~ 450 m/s, followed by another peak at f ~ 595 
m/s, and by a large peak at v ^ 800 m/s. 

Comparison of the frictional peaks of Fig. [5] with the two plateaus and jumps 
found by force pulling in Fig. [3] leads to a clear identification. The plateaus 
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are due to the rising edge of the frictional peak. So long as velocity grows, 
even marginally, for growing force, the system is mechanically stable. As soon 
as the force exceeds the peak value, however, the velocity should decrease for 
increasing force, and that causes a mechanical instability. The instability is 
then resolved by a jump to a higher stable branch. 



4 Discussion 

The main result of this work is the observation that pulling the inner nanotube 
with a well defined force of increasing (or decreasing) value generally leads to 
anomalous frictional sliding, characterized by plateaus and by jumps of the 
inner tube velocity. 

That behavior can be traced to the presence of peaks in the force-velocity 
characteristic. The frictional peaks correspond to resonant excitation of modes 
of the two-wall nanotube, when in the sliding configuration. According to 
Tangney et al. [12] , who recently found one such peak, these should be modes 
of the outer nanotube at the washboard frequency with group velocity ~ v 
(the sliding velocity). Here, however, the washboard frequencies are 3.7 THz 
and 6.3 THz for the 450 and 780 m/sec plateaus, respectively, whereas we 
observe, by analysing the velocity Fourier transforms (not shown), that the 
outer tube has large resonances at 0.14 THz and 0.65 THz, respectively. Thus 
the simple prescription seems not to work, and we must defer a detailed study 
of this aspect to further work. 

The mechanical instability under external force pulling is a novel aspect here, 
although actually quite a familiar one in a more general context. In electrical 
engineering, for example, an instability of a similar kind actually presents itself 
in all I-V systems, characterized by S-shaped negative differential resistance. 
A simple analog here would be conduction through a gas, where the current 
remains the same for a large plateau of voltages, to jump suddenly up at the 
onset of current-induced ionization, and consequent discharge. 

The phenomena presented in the present paper seem quite ubiquitous, and 
rather insensitive to many of the details of the simulations, and are therefore 
likely to occur in experiments. 
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Fig. 1. (Color online) Sketch of a F — v characteristics (typically S shaped, i.e., 
with peaks sitting on the ordinate axis), and a typical v versus F curve obtained by 
ramping F (thick solid line), with plateaus in v followed by abrupt jumps (vertical 
arrows) . 
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Fig. 3. Speed jumps and plateaus of the inner tube in a commensurate (5,5)-(10,10) 
case, when the inner tube is pulled with a constant force F = 0.12 meV/A per atom. 
Here T = 300 K. 
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Fig. 4. Chiral nanotube (11,2) inside the (12,12) nanotube pulled with a constant 
force F = 0.013 meV/A per atom. The speed increases by jumps. The motion is 
accompanied by a rotation induced by the chirality of the inner tube. 
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Fig. 5. Simulations at constrained sliding speed for the commensurate case. The 
friction force is negligible everywhere except around some velocities corresponding 
to the speed plateaus observed in the simulations at fixed pulling force. 
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